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Abstract. For a second-order elliptic equation in divergence form we investigate conditions 
on the coefficients which imply that all solutions are Lipschitz continuous or differentiable at a 
given point. We assume the coefficients have modulus of continuity satisfying the square-Dini 
condition, and obtain additional conditions that examples show are sharp. Our results extend 
those of previous authors who assume the modulus of continuity satisfies the Dini condition. 
Our method involves the study of asymptotic properties of solutions to a dynamical system 
that is derived from the coefficients of the elliptic equation. 
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0. Introduction 

We consider the regularity of weak solutions of a linear uniformly elliptic equation in divergence 
form in an open set U of R" for n > 2: 

(1) Cu :^ di{aij{x)dju) = in f/, 

where we have used the summation convention, the bounded, measurable, real- 

valued functions, and by a weak solution of (1) we mean that u € i.e. Vm is locally 

square-integrable, and satisfies 



(2) / fly {x) dj ud^r^dx^Q for aU r^eC^iJJ). 



u 

The classical results of De Giorgi [5] and Nash [20] show that u is locally Holder continuous in U . 
When the coefficients are continuous in t/, then it is well-known (cf.pj) that Vm € Li^ioJP) f*^'^ 
1 < p < oo; in fact, this is even true when the coefficients are in VMO (cf.[6]). If the coefficients 
are Dini-continuous in [/, then u is known to be continuously differentiable (cf.[13],[23]). In the 
present paper, we find conditions on the coefficients a^-, milder than Dini-continuity, under which 
u must be Lipschitz continuous, or even differentiable, at a given point. 

Let us fix an interior point of J7, which for convenience we shall assume is the origin, x = 0. 
Using a change of independent variables, we may assume that aij(O) = 5ij. Suppose that 

(3) sup \aij(x) — 5ij \ < uj{r) as r 0, 

\x\—r 

where a;(r) is a continuous, nondecreasing function for < r < 1 satisfying aj(0) = 0. We shall 
not require the Dini condition on w, i.e. r~^uj{r) £ L^(0,1); instead we assume that uj satisfies 
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the square-Dini condition: 

(4) f Lo^ir) -<oo. 

JQ ^ 

However, examples show that additional conditions are required to ensure that a solution is even 
Lipschitz continuous. 

Our additional conditions for regularity are derived from a dynamical system that we shall 
now describe. Let 

(5) R{r) / {A{re) - nA{re)e ^ 9) dse, 

where the slashed integral denotes mean value, A = (oij), r = \x\, 6 — x/\x\ e S*""^, A9 ^ 6 is 
the outer product of the vectors AO and 6, and ds denotes standard surface measure on S*""^. 
Note that < caj(r), where we use | • | to denote the matrix norm. Also note that R need 

not be symmetric. Let us consider the dynamical system 

(6) ^+i?(e-*)0 = O forr<i<oo, 

where t — ~ log r and T is sufficiently large. We shall find that the regularity of weak solutions 
of (H]) is determined by the asymptotic behavior as t — > oo of solutions of ([6]). We say that 
([6|) is uniformly stable as t — ^ oo if for every e > there exists a 5 = (5(e) > such that any 
solution of ([6]) satisfying < 5 for some ii > satisfies \4>{t)\ < £ for all t > ti (cf. 

[4]). In addition, we are interested in the condition that every solution of (0j is asymptotically 
constant, i.e. (j){t) — > (poo as t — > oo. These two stability conditions are independent of each other 
(cf. Section 5). On the other hand, it is easy to see that r^^Rir) G L^(0,e) implies that ([6]) 
is uniformly stable and every solution is asymptotically constant as t ^ oo\ in particular, if ui 
satisfies the Dini condition, then these conditions are met. 

We are now in a position to state the main result of this paper; since we are only concerned 
with regularity at a; = 0, the coefficients are not required to be continuous elsewhere. 

Theorem 1. Suppose that Oij satisfy fj) where lj satisfies and that is uniformly stable. 
Then every weak solution u £ H^^^{U) of (QJ) is Lipschitz continuous at x = and satisfies 

/ \ 1/2 

(7) \uix)-u{0)\<^IJ^^ ^^Hy)\'dy] for\x\<r/2, 

where r is sufficiently small. In addition, if every solution of is asymptotically constant, then 
u is differentiable at x — and we have 

(8) dju{0) = lini u{r9) 9j dsg. 

In this theorem and throughout this paper, c denotes a constant whose value may change with 
the instance but does not depend upon the solution u or the parameter r. 

Remark 1. // the Oij are radial functions, then R{r) = and we only require and 01) to 
conclude that weak solutions are differentiable at x = 0. Moreover, if aij{x) = a^j{\x\) + aj~{x), 
then the R in 0j is completely determined by ajj; for example, if the a\j are Dini continuous 
then weak solutions are differentiable even though aij need only be square-Dini continuous. 

We also investigate specific analytic conditions on the coefficients that imply the desired 
asymptotic properties of (jS]). Let us introduce the symmetric matrix S = —-ki^ + i-^- 



(9) S{r):^r {^[A{r9)9®9 + 9®A{r9)9]-A{r9)^dse, 
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and 

(10) l^'iS) — largest eigenvalue of S. 

In Section 4 we use the theory of dynamical systems to show that if there exist positive constants 
e and K so that 

(11) / ' ^i {S{p)) ^ <K for all s > r2 > n > 0, 

J ri P 

then ([6]) is uniformly stable. As a consequence, Theorem 1 implies the following: 

Corollary 1. Suppose that (0), and ill]) are satisfied. Then every weak solution u of is 
Lipschitz continuous at x = 0. 

What about conditions for differentiability at a; = 0? As already observed, r~^R{r) e L^{0,e) 
is sufficient, but is there a weaker condition? Let us suppose that for r G (0,£) the improper 
integral 



(12a) 



I R{p) — converges (perhaps not absolutely). 
Jo P 



Examples show (see Section 5) that this condition is not sufficient to ensure that ([6]) is uniformly 
stable; we shall require an additional condition such as 

(12b) S r R(^p^±^L\0,e), 



P 

which is also weaker than assuming R{r) r-^ e Li(0, e). In Section 4 we show that (fT2a|) and (|12b[) 
together imply not only that ([6]) is uniformly stable but asymptotically constant. Consequently, 
Theorem 1 yields the following: 

Corollary 2. Suppose that (0) and are satisfied, as well as both U2a\) and Ill2b\) . Then every 
weak solution u of (QJ) is differentiate at x — 0. 

Remark 2. Just as fl2a]) and il2b\) have replaced the more restrictive R(r)r^^ e L^(0,£), the 
assumption II 12b]) may be replaced by 



(13a) / Rip) ( / -R(o') — ) — converges (perhaps not absolutely), 



a J p 

(13b) ^fRip)(rRia)'-^Y-^eL^{0,e). 

r Jo \Jo J P 

This process may be iterated to obtain further refinements. 

Remark 3. The condition hi 2a]) can be expressed as a volume integral (computed in the sense of 
Cauchy principal value): 

A\ f ( A, \ A(x)x X \ dx . , 

(14) / A[x) ~ — i — i — ® 1 — r 1 — ;— converges for r G (0, e). 

J\x\<r \ m \x\J 

This form of the condition is better suited for changes of coordinates, so can be expressed without 
the simplifying assumption aij{0) = (5^; however, HI 2b]) is not so easily handled in this way. In a 
similar spirit, the following condition 

is sufficient for Corollary [H and easily generalizes to the case aij (0) ^ Stj ; however, it implies 



r 



-1 



R{r) £ L^{0,e), so is less general than assuming 11 12 a]) and ^12b]) . 
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Let us consider one more consequence of Theorem 1. In Section 4 we show that 

(16) [ fi{S if))— -oo asr^O 

Jr P 

imphes that the nuh sohition of ([6]) is asymptotically stable. Thus Theorem 1 yields the following: 

Corollary 3. Suppose that (0), 0), and ill}) are satisfied. Moreover, assume I116\) . Then every 
weak solution u of (QP is differentiable at x — and all derivatives are zero: dju{0) = for 
j = l,...,n. 

In Section 5 we discuss an example which illustrates the sharpness of the conditions in Theorem 
1 and its corollaries. 

Now we describe the main ideas of the proof of Theorem 1, which is given in Section 3. We 
write u in the form 

(17) u{x) = Uq{\x\) + v{\x\) ■ X + w{x), 
where 

(18) uoir):=f u{re)dsg, Vk{r):=-f u{re)Bkdse, 
and w has zero spherical mean and first spherical moments: 

(19) / w{re) dse^O^r w{re)0i dse for i = 1, . . . , n. 

We shall find that v satisfies a second-order differential system depending upon uq and w, but it 
is equivalent to a first-order system that only depends on w. Moreover, in this first-order system, 
the behavior of both v{r) and rv'{r) is controlled by the asymptotic properties of the solutions 
to To be more specific, we need to assume that a;(r) does not vanish faster than r as r — > 0: 

(20) aj(r) r"^"*"" is nonincreasing for r near and some k > 0. 

Then the assumption that © is uniformly stable ensures not only that v{r) and rv'(r) are 
bounded as r — >■ 0, but that |uo('') — '«o(0)| and \w{x)\ are both bounded by coj{r) r, so we obtain 

(21) u{x) ^ u{Q) + v{\x\) ■ X + 0{uj[r)r) as r = |a;| -> 0, 

which confirms that u is Lipschitz at a; 0. If we also know that all solutions of ([B]) are 
asymptotically constant, then v{r) — v{Q) -f o(l) as r 0, which shows that u is differentiable 
at a; = 0. 

We observe that the square-Dini condition has been encountered by several other authors 
in a variety of contexts. It was used by Stein and Zygmund [22j in their investigation of the 
differentiability of functions, by Fabes, Sroka, and Widman [3] in their study of Littlewood- 
Paley estimates for parabolic equations, and more recently by several authors (cf. [8], [3], [lO] . 
|14p investigating the absolute continuity of elliptic measure and boundary conditions for the 
Dirichlet problem. We also used the square-Dini condition in [TB] and [T7] to study equations in 
nondivergence and double divergence form. In addition, we should note that the techniques used 
in this paper are related to, but independent of, the asymptotic theory developed in [15] . 

Finally, we mention that the techniques and results of this paper apply to weak solutions of 
more general linear equations than (jlj : lower-order terms in u (even with mild singularities in the 
coefficients) as well as a nonhomogenous right-hand side (with certain integrability conditions) 
can be treated. However, sharp conditions on singular coefficients in lower-order terms requires 
additional analysis beyond the results of this paper. 
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1. Potential Theory Estimates in R" 
We will encounter an equation in the following form 
(22) - Aw^ g in M"\{0}, 

where is a distribution and we consider the derivatives in (j22p in the distributional sense. 



We will encounter certain orthogonality conditions with respect to the spherical mean, so let us 
summarize these in the following: 

Lemma 1. (a) Suppose that f e Ll^cO^"-\{0}) and g is bounded with compact support. Then 
(23a) / J{\x\)g{x)dx^ I f {x)g{\x\) dx. 



(h) Suppose that V/ G L^„^(R"\{0}). For all r > 

(23b) 7(r) = / fir9) dsg ^ ^ f 0, dJirO) dse = 0. 

(c) Suppose that V/ e L\^^{W^\{{)}) . For all r > and any i = 1, . . . , n 

(23c) / ej{re)dsg^{) ^ f dJ{re)dse^Q^ P 9,9^8, f{r0) dsg. 

J S"-i J S"-i J 

Proof. The proof of (a) is trivial. To prove (b), we consider G C^(0,oo) and compute 

1 

IS"" 







tdjdsg,(f)) ^ I /9,dif{r9)dsg(l){r)dr^ ^ / x,dj{x)(l){\x\)\x\ " dx 



^ ^ n/(x)</)(|x|)|xr"dx--^ / a;J(a;)[0(r)r^"]l=|,|0,dx 



f{r9)dsg] 4)'{r)dr = 



'0 V S"-! / 

where ' denoted d/dr. To prove (c), we again consider (f) e C^{Q, oo) and compute 



d,f{r9)dscj,{r)dr = -^ / aj(a:)|x|i-"^(|x|) 



IS"" 



y(x)[/-"0(r)]l=|,|0,dx 
f{r9)9, ds [r^-''(l){r)]'r''-^dr = 0. 

/o .^S"-i 

The proof of the remaining identity in (|23bl) is similar. □ 

Note that (I23a|) enables us to define the spherical mean of a distribution. In fact, for / g 
Li^^(R"\{0}), let us define 

(24) f{r9)^ = f(r9) - Pf{r9), 

where Pf is defined by 

Pfir9) = f fir(f>) ds^ + n0kf ^kfir(f>) ds^. 

Using 

9k9( dsg = - 5m for fc, £ = 1, . . . , n. 
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it is clear that (for each r > 0) P is the projection of / onto the functions on 5" ^ spanned by 
l,6i, . . . ,9n- For the same reason as ()23a|) , we have 



(25) / P{f)9dx= f fP{g)dx. 

for g bounded with compact support. This also allows us to define P on distributions. In fact, 
one particular instance of (|22|) that we are interested in is 



(26) -^w = [d,h]^ in M"\{0}, 

where fi G LJq^(]R"\{0}). We will solve (pS)) by convolution with r(|a;|), the fundamental solution 
for —A, but we are interested in controlling the growth of the solution near x = Q using mean 
values over annuli. We consider L^-means for any p € (1, oo): 

i/p 

(27) Mj,{w,r)^[y \w{x)\^dx 



where Aj. = {x : r < \x\ < 2r} and w may be scalar or vector- valued. To control the growth of 
the first derivatives of functions we introduce 

Afi,p(w, r) = rMp{Ww, r) + Mp{w, r). 
Proposition 1. Suppose n>2,p£ (l,oo), and f = (fi) with fi E L^oc(^^\{^}) satisfies 

|x| |/(a;)| < oo, and / |/(x)| |a;|~"~"'" da; < cx). 

|a:|<l "'|a:|>l 

Then convolution by T defines a solution w G H^^^{M"\{0}) of \26]) that satisfies Pw = and 

/ /T poo 

(28) Mi^p{w,r) <c(r'^ Mp{f , p) dp + j Mp{f,p)p-^dp 

In (|28p and throughout this paper, c denotes a constant; in other instances, the value of c may 
change line by line without change in notation. 

Proof of Proposition 1. We may assume that fi G Co(K"\{0}) since the general case may be 
handled by an approximation argument. Using (j25p . let us write the solution of (j26l) as 



w{x) = / r(|x - y\) (dd^iy) - P{dJ,){y)) dy 

[r{\x - y\) - P{r,){\y\,y)] d,f,{y) dy, y = y/\y\, 



where Tx{y) = T{\x — y\); clearly Pw — 0. To calculate ^(r^,), we use an expansion of T in 
spherical harmonics. Let Hk denote the spherical harmonics of degree k and let N{k) = dim^Hfe. 
For each k, choose an orthonormal basis {(j)k,m : m = 1, . . . ,N{k)} for Hk', for fc = 1, note that 
0i,m(^) = V^Sra- FoT uotatioual convenience, let us assume n > 3; the case n = 2 is analogous. 
For I a; I < |?;| we can write r(|a; — y|) as a convergent series 

oo I .fc N(k) 

(29) T{\x -y\)^Yl I n-2+k E (^) ■ 

A:=0 m=l 

Since / 4'k,m{0) dsg = for fc > 0, the spherical mean of F is given by 

ri\x-y\)dsy = Cn\y\^-" = Ti\y\). 
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We can also use (|29l) and the orthogonality of the 4>k,t to compute 



s"-i \y\ 



Consequently, 



By symmetry, it is clear how to modify these projections for |a;| > |?/|, so we obtain the following: 



w{x) 



\v\>\A 
f 



T{\x-y\)~T{\y\)~n^^T{\y\)x-y 



ri\x-y\)-T{\x\)-n^-^ri\x\)x-y 



d_ 

\v\>\x\ dyt 
_d_ 

\v\<\x\ % 



Ti\x-y\)-Ti\y\)~n^-^^Ti\y\)x-y 
T{\x-y\)-r{\x\)-n^-^r{\x\)x-y 



diMy)dy 
diMy)dy 
fi{y)dy 
fi{y)dy, 



where we have used the divergence theorem (and the fact that fi is supported in E"\{0}). 

If we assume r < \x\ < 2r and introduce the annulus Ar — {x : r/2 < \x\ < 4r}, then we can 
split up the integrals as follows: 



w{x) 



7^r(|x - y\)My) dy+ f A fr{\x\) + r(|a;|) x-y] Mv) dy 

A,- CfXi Jr/2<\v\<\x\ oyi \ \X\ 



d 



T{\y\)+n[]T{\y\)x-y)My)dy 

\x\<\v\<4r oyi V \y\ J 



d_ 

y\<r/2 dyi 

_d_ 

1/1 >4r (^Vi 



r{\x-y\)-n\x\)~n\^T{\x\)x-y 



T{\x~y\)-T{\y\)-n^-^^T{\y\)x-y 



My) dy 
My)dy. 



Using ([29|) we can estimate the last two integrals: 



_d_ 

\y\<r/2 dyi 



r{\x-y\)-T{\x\)-n^-l^r{\x\)x-y 



My) dy 

< cr-" 



<cl \xn\yf{y)\dy 

\y\<r/2 



\y\\fiy)\dy 



y\<\x\ 



d_ 

y\>ir (^Vi . 



T{\x-y\)-n\y\)-n'-^^V{\y\)x-y 



My)dy 



<cl \x\'\yn''\f{y)\dy 

\y\>ir 



y\>\x\ 



\yr-V{y)\dy. 
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We can also easily estimate 

lr/2<\y\<\x\ (^Vi 



T{\x\) + rMT{\x\)x-y\f,{y)dy 



< c 



r/2<\v\<\x\ 



\x\'-"\f{y)\dy 



\y\\fiy)\dy<cr-^^ \y\\fiy)\dy 



r/2<\y\<\x\ 



and 



^iri\y\)+nf^r{\y\)x-y]My)dy 



x\<\y \ <4r 



^y^ 



\y\ 



< cr^ 



< cr' 



\yn-'\fiy)\dy 

|2;|<|iy|<4r 

\yr''\f{y)\dy. 

x\<\y\ 



We conclude that 

w{x) 

Similarly, we can show that 



< c r 



/ ^T{\x-y\)My)dy 



\y\<M 



\y\\f{y)\dy + r^ / \yr^-'\f{y)\dy 



y\>\x\ 



dxjdxi 



■T{\x-y\)My)dy 



is bounded by 



y\<\x\ 



\y\\f{y)\dy + r \yr'-V{y)\dy 



Now we use Stein's inequality [21 to conclude 

/ ^Ti\x-y\)My)dy 



Lv{Ar 



< C?'|l/|lLP(i^) 



and the L^-boundedness of singular integral operators to conclude 

92 



lAr dxjdxi 
Putting this all together, we obtain 



T{\x~y\)My)dy 



<c\\f\\ 



L'-iA,.) 



Lv{Ar)- 



(30) A/i,p(«;,r)<c rMp(/,r)+r-" / \y\\f{y)\dy + \I{y)\\yr'~' dy 



y\<r 



\y\>r 



where Mp denotes the mean value over Ar instead of Ar. 

The integrals in pop can be estimated in terms of Mp and combined with the Mp term, 
fact it is elementary (cf. [16]) to establish 



y\<r 



\y\\f[y)\dy<c Mp{f,p)p"dp 



I \y\>r 

In addition, it is easy to see that 

nir 

Mp{f,r)<c Mp{f,p)p-^dp<c 

Jr/i 



\f{y)\\y\-''-'dy<c / Mp{f, p)p-^ dp. 



/4 



Mpif,p)p''dp + r Mp{f,p)p-^dp 
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Using these inequalities, it is clear that (pO| implies (|28|). □ 
Another instance of ([22|) that we are interested in is 

(31) -Aw= [f]^ in M"\{0}, 

where / £ Ll^^{M.'"'\{0}) . To control the growth of the second derivatives of functions we use 

M2,p{w, r) = r^Mp{D^w, r) + Mi^p{w, r), 

where D^w denotes the Hessian matrix of all second-order derivatives of w. The proof of the 
following is analogous to that of Proposition 1. 

Proposition 2. Suppose n>2,p£ {l,oo), and f G -^£oc(^"\{0}) satisfies 

\x\'^ \f{x)\dx < oo, and / |/(a;)| < cx). 

|a:|<l J\x\>l 

Then convolution by T defines a solution w G i?^^^(R"\{0}) of \31]) that satisfies 

(32) M2,p{w, r) < c (^r-" Mp{f, p) dp + Mp{f, p) p-^ df^ . 



2. Stability Properties of Dynamical Systems 

Let us now consider the result that we require from perturbation theory for systems of ODEs 
onT <t < oo. Without loss of generality, we assume T = 0. First, let us introduce a positive, 
nonincreasing continuous function e{t) satisfying 

/■oo 

(33) / e^{t)dt < oo. 

Jo 

Now consider the 2n x 2n system on (0, oo) 

i{t)-{: -L)e)-<')Ct)^^(o. 

where i) 7^ is a 2n x 2n matrix of the form 

(34b) n{t)^(^fj^^ f^^^ with|i?,(t)| <e(t)onO<t<oo, 

and ii) g — (.gi, 52) with gi G ^^(0, 00) and there exists (5 > so that for any choice of a G [n — S, n) 
there is a constant Cq, so that 



/oo 
\g2{s)\e-°"' ds < c„ e{t) ioi < t < 00. 



(With regard to convergence at infinity, p4cp is weaker than assuming 52 G L^(0,oo).) In 
addition, we assume asymptotic conditions on the solutions of 

(35) + = for t > 0, 

dt 

and that tp satisfies the "finite-energy condition" 

/•oo 

(36) / (|VP + |Vtne""*rfi<oo. 
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Proposition 3. Suppose that i35\) is uniformly stable. Then all solutions {(p, ^) of {S^^ that 
satisfy ^3(j\) remain hounded as t — > oo, and ijj(t) — > 0. In fact, for a ^ n ^ S with d > 
sufficiently small, we have the estimates 

(37a) sup |(/.(t)| <c(c„ + 10(0)1 + ||.gi||i), 

0<t<oo 

and 

(37b) |V'(OI <c£(i)(c„+ sup |0(t)|). 

i<r <oo 

In addition, if all solutions of i35\) are asymptotically constant as f — > oo, then the solution (0, ip) 
of ^34^ also has a limit: 

(38) (0(t),iA(i)) ^ (000,0) ast^(x,. 



Proof. Let us simplify notation by denoting d/dt by dot: d(p/dt — 0. Therefore, we want to 
study solutions of 

(j> + Ri(j) + R21P = g 

(39) 

ip — Wip + i?30 + Riip — h, 

when g £ i^(0,oo) and h satisfies the condition on 52 in p4cp for a certain value of a that will 
be specified below. Let $ denote the fundamental matrix for + i?i</i = on < > 0, i.e. 

The assumption that ([55)1 is uniformly stable is equivalent (cf. 4 ) to 

(40) \<l>(t)<l>-\s)\ < K tor t > s > 0, 

where if is a constant. Next let 5* denote the fundamental matrix for tp + R/^ip = on t > 0, i.e. 

^r + i?45f = 0, *(0) = /. 
Since £{t) — >■ as i — >■ 00, for fixed < (5 < 1 we can find ti so that 

/•oo 

(41) £(t) <5 iort>ti, and / £^{t) dt < S. 

Jti 

Without loss of generality, we can assume ti = 0. Using Gronwall's inequality, we can show 

(42) |*(t)^'-i(s)| < e'^l*-'*! for t,s>0. 

However, we also need a lower bound on ^'(t). To derive this, let tp{t) = '^{t)ipo andp(t) = |-i/;(t)p. 
Then 

p =2ipiifi H h 2lp„li>ri 

= - 2riixjjf - 2ri2'ipiip2 2ri„Vi'0n 

- 2r2iip2'(pi - 2r22ipl 2r2n'ip2lpn 



2r„i-0nV'l Srnnf/'^ 



> - 2nS\xp\^ = -2nSp. 
Integration yields p{t) > po e^^""^*; in other words 
(43) l*(OI > e~"**. 
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Now we can use the "variation of parameters" formula to conclude from the first equation in 
(1391) that 



(44) 



and from the second equation in (j39]) that 



(45) 



^-(0)+ / *-1(t)[/i(t) -i?3(r)0(r)]e-""dT 
In order to have p6p . we see from that we must have 

/•OO 

m = - / W[/i(t) - i?3(r)</)(r)]e'"^ dr, 

Jq 

and consequently (|45p can be rewritten as 

/OO 
'f-\T)[R3{T)(f>{T) - /i(T)]e-"^ dr. 

If we plug into pi)) . we obtain 



(47a) 
where 



cj^it) + s^{t) = at) = Ut) + 6 it) + 6 {t), 



pt poo 

S4>{t) ^ -^(t) j $-i(T)i?2(T)e"^*(T)y *-i(CT)i?3(cr)(^(a)e-"'"dcrdr, 



(47b) ^o(t) = Ut) = '5(0 / <^'\r)g{T)dT 



6(0 = *W/ $"^(r)i?2(T)e"^*(T) / *-i(cr)/i(cr)e-"'^d(7dr 



We want to use (|Tf|) to conclude that is bounded. 

Let X = C[0,cx)) with ||0||x supo<f<oo l'/'(i)l < oo. Notice that (gOl) and g G i^OjOo) 
imply that e X. To show ^ e X, let us use (gOl), (02]), and (p4c)) : 

pt poo 

Mt)\<Kj £(T)e("-*)"y |/i(a)|e(*-")"d(7dr 

pt pOO pt 

<K £(r) / \h{a)\dcrdT < Kca / e^ir) dr < K c^S, 
Jo Jt Jo 

although we do not care if this is small. Now let us show that S : X —i' X with \\S\\ < 1 if (5 is 

small. In fact, assume \\(j>\\x £ 1- Then 

/^t />oo 

\S(t){t)\<K I £(T)e("-*)" / e(*'")"e(a)|0(a)|dadT 

K5 



pt pOO pOO 

<K / £(T)e("-*)" / e(*-")"£(a)dadT < / e^{T)dT< 

Jo Jt n-S Jf^ 

so II 5*11 < 1 if (5 is sufficiently small. We conclude that = (1 + 5)^^^ G X, i.e. (j> is bounded. 
Now if we apply (gH) to (|46]), we find that 



1^^(01 < ce("-*)* / e(*-")"£(r)|0(T)|dr + e("-^)* / |/i(T)|e(^-")" dr 



< c£(i) ( Cc + sup |0(t)| ) , 

i<r<cxD 
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which is ()37bp . We can then combine this with (j44|) and (|41|) to conclude 

m)\ < c {\m\ + Ml) +cs{c^ + mx). 

Taking 5 sufficiently small, we can conclude (|37ap . 

If all solutions of (1351) are asymptotically constant as t — >• oo, then $(oo) — limt^oo ^(t) exists 
and from (l44l) we find 



0OO = H^) + J $-^(r)[g(T) - R2{T)4,{T)]dT ^ 

To show (j){t) (/)oo we estimate three terms: |($(t) - $(cx))) 0(O)| < |$(i) - $(oo)||(/)(0)|, 

$(0 / $"l(T)c/(T)dT - $(CX)) / <^-^{T)g{T)dT 

Jo Jo 



< c (^\m-H^)\\\9\\i+ \9iT)\dT 



pt POO 

/ (T)i?2 (T)V'(T)dr - $(oo) / $-i(T)i?2(r)?/'(r)dr 

<c{l + Mx) 



/OO 
e2(T) dr 



This confirms dSSl 



□ 

It will also be useful to know that the uniform stability is not affected by perturbation with 
terms bounded by e^(t). 

Proposition 4. If R{t) — R{t) G L^{0,oo), then ^ is uniformly stable if and only if the same 
is true of 

d(j} 



dt 



+ R<j> = 0. 



The same equivalency applies if the property that all solutions are asymptotically constant is added 
to uniform stability. 

Proof. Let S — R — R, denote the fundamental solution for ^ + Rcj) = on t > 0, and 
$ denote the same for ^ + R(p = 0. According to we know that (PO]) holds and we want to 
prove that a similar bound holds for But we can solve 



dt 

by variation of parameters to obtain 



+ R(j>= -Scj), 



(f){t) = $(i)$-i(s)e - / $-1(ct)S'(ct)0(o-) da, 



Applying Gronwall's lemma and (|40l) . we obtain 



\(j){t)\ < K eiip (^k \Sia)\da^\^\< K exp(^k \Sia)\da^ \^\. 

But (j){t) = <I>(t)$~^(s)^, so this last estimate shows ^ + R(j> = is uniformly stable. The 
additional property that solutions are asymptotically constant also follows from the variation of 
parameters formula. □ 
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3. Proof of Theorem 1 

Recall that aj(r) is a continuous, nondecreasing function for < r < 1 satisfying uj{0) — 0, 
and (I20p . Since our result is local, we focus our attention on a very small ball centered at a; = 0; 
by rescaling, we may assume this is Bi{0), the unit ball centered at a: = 0, and that ([3]) holds for 
< r < 1; in fact, given any small 6 > 0, we can similarly assume that 

(48) < S and uj{l) < S. 



If the Gij are continuous in U, then we know by the results of that Vu £ L^od^) ^"-"^ 
p G (1, oo). More generally, under the small oscillation assumption (jS]), we can fix any p € (1, c») 
and conclude that Vu G L^^d^) provided 5 = 5{p) in is sufficiently small; cf. Corollary 
6.2 in [18]. Henceforth, we fix p > n and assume that 5{p) has been chosen small enough that 
Vw G LP(Si(0)). 

For our analysis, it is advantageous to extend the problem to all of M", so let us redefine and 
extend the aij outside of i?i(0) by 

(49) <^ij{x) = ^ij for I a; I > 1. 

If we extend cj(r) to be u>{l) for r > 1, we see that Lu{r) is still nondecreasing and uj{r) r~^^^ is 
still nonincreasing; in particular, we easily see that 

(50) Mpiujf, r) <cuj{r) Mp{f, r) for < r < cx3, 

where / G Hg^^{W"-\{0}) and c is independent of r. Now let us introduce a smooth cutoff function 
x{r) which is 1 for < ?' < 1/4 and for r > 1/2. Given a weak solution u of ([1]), we can write 

di{aijdj{xu)) = di{aijx'Oju) + x'f^iO-ijdjU. 

Since we are interested in the behavior near a; = where xd^^D'^^C^^) ^-nd u{x) agree, we can 
assume that u is supported on |x| < 1/2 and satisfies 

(51a) d^{a,j{x)dju) = dj^ + fa in K", 

where /i, /o G L''(M") are both supported in 1/4 < \x\ < 1/2, and (using ([2]) with rj — x) 



(51b) / fo{x)dx = Q- 

when convenient, we let / = (/i, . . . , /„). Of course, we now must replace ^ with the following: 
(52) / aij{x)djudi'ridx = / fidirjdx- I foridx, 



for all rj G C^(K"). 

Recall the decomposition u{x) = ua{r) +v{r) -x + wlx) given in (|17p . Orthogonality properties 
show that 



(53) Wue L^{Bi{0))^ f {{u'^f + \v\^ +r^\v'\'^)r"-'dr <oo and Wwe L'^{Bi{0)). 

Jo 

We want to show that |uo('') — wo(0)| ^ and |u'(x)| |a;|^"'^ as r = |a;| 0, so that the 
differentiability of u is determined by the behavior of u as r — )■ 0. The strategy is to show that 
V satisfies a system of ODEs which depends upon w and that w satisfies a PDE which depends 
upon V. We obtain the system of ODEs by plugging u{x) ~ uo{r) + v{r) ■ x + w{x) into ([52]) and 
choosing special rj G C^(R"). For example, taking rj — r]{r) G C^[0,oo), we obtain 



(54a) / (^a + r/3 • u ' + 7 • u + p[Vw]^ rj' r" ^ dr = j (^f 1]' - fovj r 



dr. 
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where 



J 5"-i J S"-i 



(54b) lj{r) = r a,j{r0)0,dsg, p[Vw]{r) ^ / a,j{r6)djw{r6)e,dse, 



/.(r-0)^.dse, k{r)^r k{re)dse. 

Using ([3]) and Lemma 1, we see that these terms satisfy 

\a{r) - 1|, \l3{r)l \^{r)\ < w(r) for < r < 1, 



(54c) |p[Vw](r)| < uj{r) f | Vw| ds for < r < 1, 

/('■) =^(0 - for < r < 1/4, 

whereas for r > 1 we have a{r) — 1 and all the other terms vanish. We can integrate (|54a|) to 
find 

(55a) Q;(r)uo(r) + r/?(r) • v'{r) + j{r) ■ v{r) + p[Vw](r) = i^(?'), 

where 

(55b) z?(r) = /(r) + r^-" /' J,{p) dp 



has support inl/4<r'<l/2. (Note that i?(r) = for r < 1/4 follows from (I54ap with supp 
?7 C [0, 1/4) and r?(0) = 1, whereas i){r) = for r > 1/2 is a consequence of (I51bp .') 

Similarly, we can let 77 — r]{r)xi in ([52]) and obtain a 2nd-order linear system of ODEs. We can 
use ()55ap to eliminate mq and then reduce the 2nd-order system for t; to a Ist-order system for 
{v, Vr)] of course, these systems also depend on w. This Ist-order system is simplified by changing 
independent variables io t = — logr, so we introduce 

(56) e{t)=u:{e-'), 

which satisfies ([33l) by (|4]). In the Appendix, we show that the Ist-order system for {v,vt) may 
be converted to the form (|34p : 

s(v*) + (S ( J) - 'ft '-' + '■(0, 

where 7?. = for i < 0, but for t > it is of the form (|34bp with 

(57b) i?i(0~/ {A - nAe(^e)dse as i ^ 00, 

where ss indicates that the difference is bounded by ce^{t); the term g(t, Vw) = for t < but 
satisfies 



(57c) \g{t,Vw)\<ce{t)/ \Vw\ds for t > 0, 

J S"-i 

and the term h is in with support in log 2 < t < 2 log 2 with L^-norm satisfying 

(57d) ||/^||i<c (||/||p + ||/o||p). 

Moreover, the new dependent variables {4>,ip) are related to {v,vt) according to: 

(57e) (/^] - ^ < ceit) (\v{t)\ + WM + f\Vw\ds 



DIFFERENTIABILITY OF SOLUTIONS 
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Now, given w with suitable properties, we solve ()57a|) with initial conditions 0(0) = = ■0(0) 
to find {(f), ^) and hence v. But we want to separately control the dependence of v upon w, so let 
us write v = + v'^ where corresponds to solving (j57ap with h{t) = and v'^ corresponds 
to solving (|57ap with g{t, Viy) = 0. We want to apply Proposition |3] to estimate i/™ on (0, oo), so 
we need to confirm that g — (51,(72) satisfies gi G L^{0,oo) and g2 satisfies (|34cp . To show that 
gi £ L^, we use (|57cp to conclude 

/>oo / />oo \ 1/2 / /.oo /" \ 1/2 

/ \giit,\7w)\dt < c ( s^{t)dt] ( / iVwpdsdij 

1/2 . ^1 ^ , . 1/2 



and then invoke (|i8l) to conclude 

1 /2 

(58) \\g^\\^<cV6 (^J^y\Ww\'ds^^ . 

Similarly, to verify (|34cp , we estimate 

/oo nOO n 

\g2{T)\e"^'' dT <ce{t) e""^^-^^ w\ ds dr < c^eit) , 

where 

(59b) = I v.p ,,;;£)"'. 

As we shall see below, M2(V?i;,r) — 0{Lu{r)) as r — > 0, so the finiteness of Cq, and the bound in 
(|58p follow from the following calculations: for j = 0, 1, . . . , let = 2"-' , so that 

(60a) / / \yw\^ds^^y / \\/w\''ds^<cyM^{Vw,r,) 

Jo Js^ P j-'Jr, P 

and 

00 00 

)) Tu;Hr,) = yu;Hr,)^-^ _ , „^ 

rj+i Jq P 



(60b) 5:a;^(r,) = yu;\r, f-l^^ < c / c.2(p)^ < c5, 

where we have used (^5]) at the end. 

Now let us derive the PDE for w. Introduce the matrix 51 — (ftij) with entries 

(61) ilij = aij — Sij, 

and recall that |f7y(r)| < aj(r) for < r < 1 and ilij{r) — for r > 1. We can apply I — P to 
(|51ap . to obtain 

[a, (a,, (rr)5, u)] ^ = [d, /, + /o] ^ in R" . 

If we substitute p7|) into this and use P[A(uo + 1/ • if)] — A(mo + • x) and P[At(j] = 0, we obtain 
the following equation for w: 

Aw + [div{nVw)]^ + [div{nV{v-x))]^ + [div(17Vuo)]^ = [ft/, + fo]^. 

But we can use (|55aP to eliminate uq and v = v"' + to write this as: 

Aw+[div(17Vt«)]-L + [div(17V(i7"' • x))]^ - [d\^{a-^m{rP ■ {v"")' +^-v'" +p[Vw]))]^ 

(62) 

=[ft/* + /o + div(a-^ 176'(r/3 • (w°)' + 7 • t;") - f7V(t7° • f)) - Aw{a-^ned)]^ . 



16 



VLADIMIR MAZ'YA AND ROBERT MCOWEN 



Our strategy now is to simultaneously solve (|57ap and ()62p for v and w, then we can plug these 
into ()55ap and integrate to find uo with uo{r) — for r > 1/2. Indeed, for our chosen p > n, we 
assume w is in the space Y consisting of functions in H^^^{M."\{0}) with finite norm 

(63) \H\y ^ sup ^- +sup — - — . 

o<r<i uj(r)r r>i or " 

We plug w into (|57ap and solve as described above to find v — v'^ + v'^; since M2(Vz/;,r) < 
cMpC^w, r) = 0{uj{r)) for < r < 1, we can apply Proposition 3 to estimate w™. Now we want 
w to satisfy so let us apply A^^, i.e. convolution by the fundamental solution, to both sides 
of (El to obtain 



(64) w + T[w] = t 

where 

T[w] = ([div(l]Vw)]^ + [div(l]V(t/"' • x))]^ - [div {a' ^ ne {r f3 ■ (w™)' + 7 • w"" +p[Vw])]^) 

e = A"^ ([aj, + /o - div(a-i06lz9) + div {a- ^n9{r (3 ■ (w ")' + 7 • tjO) - f7V(i/° • x))]^^ . 

We want to use the results of Sections 1 and 2 to show ^ E Y and T is a bounded operator F — > F 
with small norm. 

To show r : y — > y has small norm, consider T[y] for < 1, i.e. we assume y satisfies 

(65a) Mi^p{y, r) < u}{r) r for < r < 1, 

and 

(65b) Mi^piy, r) <6 r"" = w(r) r"" for r > 1. 

Now let us consider separately the three terms, 

Ti[2/] -A-1 ([div(r!Vy)]^) , T2[y] - A^^ ([div(l]V(«^ • f))]^) , 
r3[y] = A-i (^[dW{a-^n9{r^ ■ (w ^)' + 7 • w ^ + p[Vy]))]^) . 
First we consider Ti. Using Proposition 1 and we have 

(66) M,^p{Ti[y],r) < c (^r"" ^ Mp{nVy, p)p" dp + Mp{nVy, p)p-^ dp^ . 

But recall |i^(r)| < uj{r) for < r < 1 and ^{r) = for r > 1. Thus, for < r < 1, we can 
use (I65ap and (j65bl) with the facts that Lu{r) is nondecreasing and a;(r)r~^+'* is nonincreasing to 
obtain 



^g^^^ Mi,p(Ti[2/],r) <c (^r-"y^ ^^(p)p"dp + r^ j ^'(p)p-'dp 

< c {r-''uj^{ry'+^ + r^Suj{r)r-^+^{r-^ - 1)) < C(5a;(r)r. 

For r > 1 we simply get 

(67b) Mi^p{Ti[y],r)<cr~" [ p" dp < cS^ r'" . 



The estimates (I67ap and (j67bp show that ||Ti[?/]||y < c6, so for S sufficiently small we conclude 
that Ti : y — > y has norm less than 1/3. 



DIFFERENTIABILITY OF SOLUTIONS 
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Next we consider T2. To use Proposition 1, we need to estimate Mp{y{vy ■ x),r). But and 
r{vy)' can be expressed in terms of the solutions -0 of (|57ap with = 0, to which we can apply 
Proposition 3 to find 

(68) sup |V(w^ ■x)\<c sup(r|(v^)'| + < c sup + IV-I) < c (c„ + ||.gi||i), 

|x|<l r<l t>0 

where is given in (j59b|) and ||(7i|| is estimated as in But, using (|60ap . p > 2, (|65ap . and 



(|60b[l . we have 

pi p rl ^ °° °° 

/ > \Vy\^ds^<cY,Kq{Vy,T,)<cY,Ml{yv,r,)<cY,u\r,)<c6. 
Jo J 5-1 P ^.^0 ^.^0 j.^o 

We conclude that for < r < 1: 

(69) Mpiyiv^ ■ x),r) < c sup |V(w^ • x)\ <cVS. 

\x\<l 

We can use this in Proposition 1 and estimate as above to obtain for < r < 1: 

(70a) Mi^p{T2[y],r) < c (^?'"" ^ oj{p)VSp" dp + uj{p)V6p-^ dp^ <coj{r)VSr. 

Meanwhile, for r > 1 we know v — = 0, so by Proposition 1 

(70b) Mi^p{T2[y],r) < cr-" f uj{p) Mpi^iv^ ■ x), p)p'' dp < c6^/^r-". 



For 5 sufficiently small, we conclude from (j70ap . (j70b|) that T2 :Y ^ Y has norm less than 1/3. 

To show Ta : r -> r is smah, we need to estimate Mpia^'^fieirP ■ {vy)' + 7 • w ^ + p[Vy]), r) 
only for < r < 1 (since it vanishes for r > 1). But recalling the bounds on fl, /3, 7 and p[Vy], 
we have for < r < 1: 

Mp{a-'^n0{r0 ■ {vyy + j-vy) +p[Vy],r) < cuj^ir) {Mp{r{vyy ,r) + Mpiv^ ,r) + Mp{Vy,r)) 

where at the end we have used Proposition 3 and (|68p , similar to our derivation of (|69p . Applying 
Proposition 1, we obtain for < r < 1 

Mi,p(T3[?;],r)<c(^r-"^ VS uj\p) p'' dp + VS uj^p) p-^ dp^ < cS^^^ uj{r) 

and for r > 1 we have simply 

MijT3[y],r) <cr-^ [ VS uj{p)p'' dp < cS^/^ r-" . 



Taking 5 sufhciently small, we conclude that T3 : y — > F has norm less than 1/3. We have 
therefore shown that T = Ti + T2 + T3 : Y ^ Y has norm less than 1 . 

To show ^ GY, first note that supp / C = {x : 1/4 < |x| < 1/2} implies Mp(/, r) = for 
r < 1/8 and r > 1/2, whereas Mp{f,r) < c||/||p for 1/8 <r< 1/2. Now we separately treat 

and 
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Since both fi and /o are supported in A1/4, we can apply Proposition 1 to estimate ^1 as 

M (C forO<r<l, 

and Proposition 2 to estimate ^2 as 

A^,,(6,.)<M.,(6,.)<{'^||fjj"'!.. 

[c||/o||pr for r > 1. 

Since nOd = for |a;| = r > 1 and r < 1/8, whereas Mp{a-^nei},r) < c{\\f\\.p + ||/o||p) for 
l/8<r<l, we similarly conclude 



ci\\f\\p + \\fo\\p)r^ forO<r<l, 
c(ll/1lp + ll/o||p)r-" forr>l. 



To estimate ^4 we need to estimate Mp{\/{v^ ■ x),r). But t/° and r(v°)' can be expressed in 
terms of the solutions 0, ^ of (j57ap with g = 0, to which we can apply Proposition 3 to find 

sup |V(u° • x)\ < c\\h\\i. 

\x\<l 

Combined with ()57dp . we find Afp(V(v ° • x), r) < c (||/||p + ||/o||p): and then applying Proposition 
1 yields 

Mi.p(a, r) < c {\\f\ + ll/ollp) (^r-" £ ^{p)p'' dp + j' Lo{p) p-^ dp 

< l'c(||/1|p + ||/o||p)^(r-)r forO<r < 1, 
" lc(||/1|p + ||/o||p)r-" forr>l. 
Plugging these estimates into we find that ^ G F with 

llellr <c(||/||p + ||/o||p). 
We conclude that (|64|) admits a unique solution w satisfying 

(71) <c(||/||p + ||/o||p), 

and then we use this to find and as described above. We also know from Proposition 1 
that Pw = 0, which in particular shows that /|^|^^ w(x) dx = for every r > 0. Since p > n, we 
can apply Morrey's inequality (cf. Theorem 7.17 in fT2l) to obtain 



(72) sup |w;(a;) I <c„r/ \Vw\p dy \ . 

\x\<r \J \y\<r J 

But for fixed r G (0, 1), we can introduce rj = r 2^^ to compute 

/oo „ 
iVwl" dy = ^ 2-^"rT" / |Vu;|p dy<c sup Mp(Vw,p). 

\y\<r ^.^Q ■/rj + i<|y|<rj 0<p<l 

Recalling that ^ implies Mp{Vw,r) < cuj{r) {\\f\\p + \\fo\\p) for < r < 1, we find 

(73) sup <crc^(r)(||./||p + ||/o||p) for < r < 1. 

|a;|<r 

In particular, this implies that w is differentiable at x = with djw{0) — 0. 
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What about v and uo? Since we now know Mp(Viy, r) < ccj(r) (H/llp + ||/o||p), we obtain 

„1 „ , oo oo 

/ >|VH^d.^ < cJ2M^iVw,r,) < c(||/|l,+ ll/ollpf ^a.2(r,) < cS i\\fX + WfoWp)' ■ 
Consequently, our analysis of (|57aP yields 

(74) sup (Kr)|+r|T;'(r)|) <c(||/||p + ||/o||p). 

0<r<l 

Using (|54ap . we perform the following estimates 

\uo{r) - uom < \uo{p)\dp<c (^mp)\ + \pP{p)-v'{p)\ + \^{p)-v\+u;{p)j^\Vw\ds^ dp 

[mp)\ + \pm ■ v'ip)\ + mp) ■ v\) dp < c^m r {\\fX + wmi) 

u:{p)j^\Ww\ dsdp < c^(r) r |Vu;| ds ^ < cu{r) r (\\fX + ||/o||) 

to obtain 

(75) ko(r) - 1.0(0)1 < cujir)ri\\fX + \\fo\\p). 

But ([75)1 implies that uq is difFerentiable at r = with Mo(0) = 0. 

Thus we have found a solution of (j51ap in the form u{x) = uo{r) + v{r) ■ x + w{x) which is 
Lipschitz continuous at x = 0, but we need to verify that u coincides with the solution u of (I51ap 
that we began with. However, if we let z{x) = u{x) — u{x), we find that 2; is a weak solution of 
the homogenous equation £2: = in R", and z{x) — >■ as x — >■ 00 (since |a;| > 1 implies u{x) = 
and u{x) = w{x) = 0{\x\~") as |a;| — 00). The maximum principle shows that z = 0, i.e. u = u. 

We conclude that our solution u of ([T|) is Lipschitz continuous at x = 0. To obtain the desired 
estimate (O, we first combine (|73p . ([74| . ([75l) . and recall the definitions of / and /o to conclude 

(76) \u{x)-u{0)\<c\x\{\\yu\\L,^B') + MLHB')) for0<|x|<l/2 and B' = ^1/2(0). 
But as a solution of (IT|), u satisfies the elliptic estimate 

(77) ||Vu|Up(s') <c\\u\\l.{b*). where B* ^ ^3/4(0), 

which can be found, for example, in |lj when the coefficients are continuous; however, their proof 
extends directly to the case where the coefficients have small oscillation, which we may assume in 
the unit ball by taking 5 sufficiently small. But from [TH], u also satisfies the following estimate: 

(78) sup \u{y)\<c\\u\\L2^^B), where B = Bi(0). 

ty|<3/4 

Using these in (I75|) . we obtain 

(79) \u{x) - um < c\x\\\u\\L2^B) for < |x| < 1/2, 

which is ([7]) for r = 1. The case of general r G (0, 1) can be achieved by scaling: x = x/r and 
u{x) — u{x). Thus ([T]) is proved. 

Now let us add the hypothesis that every solution of ([6]) is asymptotically constant. Then, 
according to Proposition 3, the solution (0, ip) satisfies 0(t) — (/)oo and ip{t) — >■ as t — >■ 00. Using 
(|57e|) . we see that — )■ as t ^ 00 and consequently 

lim v{r) = n lim (/)(i) — n(j),^ and lini rv'{r) — 0. 

r— ^0 t— >oo r— 
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In particular, w(|a;|) -x is difFerentiable at x = 0. Putting this together with what we have aheady 
found about uq and w, we conclude that our weak solution u of ([T]) is differentiable at x = and 
® holds. □ 

4. Proof of Corollaries 
To prove Corollary 1, let us write ^ as 

(80) = B{t)(p for T < t < oo, 

where B{t) is bounded, but not necessarily self-adjoint. (In fact, we know B{t) ^ as t ^ oo, 
but we will not need this fact here.) Let ii{B) denote the largest eigenvalue of {B + i?*)/2, which 
satisfies (cf. % Ch.II, Sec.l) 

\l + hB\~l 

n(B) = lim J ! . 

h-i-o+ h 

Assuming that B{t) is continuous, the following inequality is proved in f4] (Ch.III, Sec. 2): 

(81) \(t){t)\ < \(f>{s)\ exp f^{B{T)) dr^ for t > s > T. 

Let us verify that (|8ip holds even though B{t) may be discontinuous. 

Let f{t) — \4'{t)\ for t > T. It is easy to see from (|80| that (f> is Lipschitz continuous, so 
Rademacher's theorem implies that (j) is differentiable almost everywhere. Consequently, almost 
everywhere / has a right-hand derivative f+ satisfying 

Uit)^ lim m±M^^_m= hm MW±MM_MM < a.e. 

/i^o+ h h^o+ h 

Since ii{B{t)) is bounded, we know that 

w{t) ■.= f{t) exp j^^i(B{T))dT^ 

is continuous and has a right-hand derivative satisfying 

w+{t) = [f+{t) - ix{B{t))f{t)\ exp ^x{B{t)) dT^ < a.e. 

We conclude that w(t) is nonincreasing for almost every t > T. But wit) is continuous, so wit) 
is nonincreasing for alH > T. We conclude that (|8T1) holds. Moreover, ([8T|) together with 



f h{B{t)) dr < K for i > s > T 

J s 



implies that ([50)) is uniformly stable. Thus we may apply Theorem 1 to obtain Corollary 1. 

To prove Corollary 2, we assume p2|) and introduce a change of dependent variable (as in [7] 
Section 11.1): 

cj,{t) = {I + S{t))m 

where 

/oo 
i?(r)dr, with i?(t) = i?(e-*). 

We find that ^ satisfies 

{I + S{t))^+RS^ = 0. 
dt 
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But S{t) — ^ as i — > OO, SO we can take T sufficiently large and conclude that (/ + S{t)) is 
invertible, and consequently ^ satisfies 

(82) ^+Q(t)e = 0, 

where by hypothesis we have Q ~ {1 + S)^^RS G L^{T, oo). As we have observed, Q G L^{T, oo) 
implies that all solutions of ([82l) are asymptotically constant, i.e. ^(t) — £^oo + o(l) a.s t oo, so 

<f>it) = (/ + SitMoo + o(l)) - Ceo + 0(1) 

and we see that all solutions of ^ are asymptotically constant. We have also observed that 
Q G L^iT, oo) implies that (l82t is uniformly stable, so the same is true of 

To prove Corollary 3, observe that ^{B{t)) dr — ^ — oo as t — > oo together with ([5T|) implies 
that all solutions of ([5D|) tend to zero as i — > oo, i.e. the null solution is asymptotically stable. 
But this implies that v{r) in (PT|) satisfies v{r) — > as r — > and Corollary 3 follows. 

5. Examples of Gilbarg-Serrin Type 
In [11] , Gilbarg and Serrin consider examples of the form 

(83) aij{x) = + g{r)9i9j, 

where g{0) ~ but vanishes slowly as r — > 0. They use such examples to show that Dini continuity 
is essential for their "extended maximum principle" to hold, but we shall use them to explore the 
conditions in Theorem 1 and its corollaries. We assume that |.g(r)| < uj{r) for r near with uj 
satisfying Q, and we can explicitly calculate the quantities introduced in the Appendix: 

1+3(0 r ^/ ^ _ A , 5W~ 



a(r) = l + g(r), /3(r) = = 7(r), A{r) = B{r) ^ ' I, C(r) = 1 + ^ 

n \ n 

Moreover, the matrix ([S]) is given by 

1 — n 

Rik{r) = g{r)Sik, 

n 

so the dynamical system ([6]) reduces to the scalar equation 

d(t> n — 1 , , 

(84) = 5(00, 

at n 

where g{t) — g{e~*). 

Now consider a weak solution u of ([ij in a domain containing x = and aij of the form (j83[) . 
According to Theorem 1, m is Lipschitz continuous at a; = provided (|84l) is uniformly stable for 
t > T with T sufficiently large; but it is easy to solve (|84|) and see that it is uniformly stable if 
and only if 

(85) / 5(t) dr <K for t > s > T. 

J s 

Moreover, ^{TZ{r)) = (1 — n~^)g{r), so pT|) agrees with (1551) and we see that Corollary 1 is sharp 
for this class of examples. On the other hand, solutions of (|M| are asymptotically constant if and 
only if the improper integral 

poo 

(86) / ^(t) dr converges to an extended real number < oo. 
Jt 

Thus Theorem 1 implies that u is differentiable at x = if both (|85l) and (1551) hold. The case 

/"OO 

(87) / g{t) dt = -oo 
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in (|86|) pertains to Corollary 3, which is sharp for this class of examples. On the other hand, 
the case that g{t) is integrable pertains to Corollary 2 and coincides with the hypothesis (|12a|) : 
however, in Corollary 2 we also require (jl2bp . since (as mentioned in the Introduction) the 
condition (|12aP alone does not imply the uniform stability of ([6]). 

In fact, this class of examples may be used to show not only this last statement, but in general 
that uniform stability is not implied by every solution being asymptotically constant: we only 
need to construct g{t) for which (|86|) holds but (f85|) fails. Moreover, if (|86|) holds because g is 
integrable on (0,oo), then this explains the need for a condition such as (jl2bp in Corollary 2; on 
the other hand, if we construct g for which ([57| holds and yet (|55|) fails, then we see that (ITTIl 
is not implied by ([T6|. so both conditions are necessary in Corollary 3. In this latter regard, let 
us observe that [4] gives an explicit example of a function g{t) satisfying (jST)) and yet (|85|) fails: 
there exist tj oo for which J^^'^^ gij) dr — >■ c» and yet J^^^^^ 9{t) dr — oo more rapidly so 
that ([87|) still holds. Now the example in ^ does not have g{t) — > as f — > oo, but the example 
can be modified to achieve this; in fact, we can even arrange g{t) = 0{t^'^/'^), which implies that 
g G L'^iT, oo) and so the aij are square-Dini continuous at a; = 0. Moreover, the example can be 
modified so that ([87|) is replaced by the condition that g is integrable on (0, oo). Thus (|85|) and 
(|86p are completely independent conditions, even under the assumption that the coefficients aij 
are square-Dini continuous at x ~ 0. 



6. Appendix 



In this appendix we provide the details behind the derivation of the dynamical system ([571 
To express this system, let us introduce n x n matrices A, B, C and vectors ^, C by 



£,t\^'w\{r) ^ / aijOiOidjwdse, [Vw] (r) = aijl 



atjirO) dsg 



idjW dse 



which satisfy for < r < 1 



-'II \B-n-'l\, \C-I\<Oj{r) 



(88) 



|C[VH(r)|, mw]{r)\<ujir] 



iVwl ds, 



while for r > 1 we use (|49|) to conclude A — n^'l = B, C = I. (Notice that the matrix A{r) 
introduced above is not the same as the matrix A{x) used in the Introduction.) Now using 
V = Ti{r)xi in ((52|) . we obtain the 2nd-order system of ODEs 



(89) 
where 



r"(M[,/3 + rAv' + Bv + e^[Vw] - f*)] + r''-' (u'^j + rBv' + Ci/+ C[Vw] + /^) = 0, 



ffir) = 



5„-i 



h{rB)e,ease and fl{r) = 



h{rB)e, dse 



are supported in 1/4 < r < 1/2. Next we can use (|54ap to solve for Uo(r') and eliminate Uq from 
I; this is the 2nd-order system of ODEs (depending upon w) that we want to analyze: 

{rAv' + Bv + ^[Vw] - / - (r/3 • w ' + 7 • t7 + p[Vw] - i?) /3 

r" (rBv' + Cv + C\^w] + p\ - (rl3 ■ v' + j ■ v + p\\/w] - i^) 7 

V / a{r) \ I 



(90) 
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If we make the substitution r = e * , then ([90|) becomes 



-Avt + Bv + i[Ww] - f* - -(-/3 • vt +^-v+p[Vw\ ~ 

a 



a 



which after some rearrangement can be written 



'Avt + Bv + ^[VwJ — ^ p 

a 



-{B- nA)vt + ^-^(7 - n/3) 
a 



+ {C-nB)v- — (7-n^) 



p[Ww] - d 



To avoid differentiating the coefficient matrices, let us convert this to a Ist-order system for 
the 2n- vector V — (Vi, V2) where Vi = v and 



(91) 



V2 := -Avt + Bv + £,[\/w] - f 



■# 



P ■ Vt ~ ^ ■ V ~ p[Vw] + d 



Notice that, as in Section 2, we have omitted arrows over the vectors Vi, V2, and V, and we 
shall now also use the dot notation for d/dt. Since A is invertible (for S sufficiently small), our 
Ist-order system can be written as 



Vi-A-^BVi+ A-W2 + 



I3-Vi~j-Vi 



A-^(3 = A- 



,1 /# Mv^]-^ 

a 



(^2) V2+{C - BA-^B)Vi + {BA-^ - n)V2 + ^^(7 -{n + A-^)p) 



+ - -{{n + A-^)P~^)=n 

a 



■ P[V^^] - ^ g 
^[Vw] /3 



7-C[VH -/^ 



Now ([5^ implies that 

and if we use the second equation in (|92p we see that 

(\V2? + \V2?) 



We can summarize this as 
(93) 



{\V\^ + \V\ 



dt < 



Notice that the terms involving Vi and V2 in ((M)) are of the form {I + D{t))V where / is the 
identity matrix and |-D(t)| < ce^(t). So we may multiply by (/ + D{t))~^ and, after some 
calculations, see that V satisfies a Ist-order system in the form 



(94a) 



^ + M{t)V ^ F[t,Vw) + Fo[t), 
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where M{t) is a 2n x 2n matrix of the form 

M{t) = Afoo + 5i(<) + 52(t), 
-/ nl 



(94b) ^^^{{1-71-^)1 {l~n)lj' 

/ I-A-^B A-^-nV 
^^^^>-\C-BA-^B + ^I BA-^-IJ' 

the Si satisfy 

|5'iW| < eW and \S2it)\ < ce^{t) for t > 0, 
the vector F(t,Ww) satisfies 

(94d) \F{t,\'w)\ < ce{t)P |Vu;|ds for t > and F(t,Vu;) = for i < 0, 

and Fo{t) has support in log 2 <t< 2 log 2 with L^-norm satisfying 
(94e) ll^^olli < c(||/||p + ||/o||). 

We can calculate the eigenvalues of M^o to be A = (n times) and X = —n {n times). The matrix 



nl nl 
I (1-n)/ 



(95a) J 
diagonalizes Moo, i-e. 

(95b) J^^MooJ = diag(0, . . . , 0, -n, . . . , -n), 

so we introduce the change of dependent variables V — >■ {4>, ip) by 



(95c) V = J 



We find that the dynamical system (|94ap now takes the form (|57ap . where the conditions (|57cp 
and (|57dp follow from (I94dp and (|94cp respectively, and TZ is of the form (|34bp with 

71 — 1 71 — 1 1 

(96) Ri = — ^A-i A-^B + C- BA-^B + -BA'^ - I. 

To simplify this expression for let us write 

A = n-^{l + A), B = n-^{l + B), and C = n'^l + C), 
where \A\, \B\, \C\ < ce{t) a.s t oo. Then 

A-^ w n{I - A), 

where w means that the difference is bounded by ce^(t) as f — cx), and a calculation shows 

RiKC-B = C-nB as i cx). 

Since 

C-nB=f {A-nAe(E)0)dsg, 

J S"-i 

we see that (j57bp holds and we have completed our derivation of (f57|) . □ 
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